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The understanding of the interaction of nucleons in nuclear and neutron-rich matter at non-zero 
temperature is important for a variety of apphcations ranging from heavy-ion coUisions to nuclear 
astrophysics. In this paper we apply the Dirac-Brueckner-Hartree-Fock method along with the 
Bonn B potential to predict single-particle properties in symmetric nuclear matter and neutron-rich 
matter at finite temperature. It is found that temperature effects are generally small but can be 
significant at low density and momentum. 
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I. INTRODUCTION 

The nuclear equation of state (EoS) at finite temper- 
ature is of fundamental importance for heavy-ion (HI) 
physics [l| as well as nuclear astrophysics, particularly in 
the final stages of a supernova evolution. Knowledge of 
the finite-temperature EoS can be of great help in the in- 
terpretation of experiments aimed at identifying a liquid- 
gas phase transition. Presently, findings concerning such 
transition are very model dependent [^tSI- 

When other aspects, such as spin asymmetries of nu- 
clear/neutron matter are included as well, conclusions 
are even more contradictory. For instance, the influence 
of finite temperature on the manifestation of ferromag- 
netic instabilities is an unsettled question. In Ref. 
phenomenological Skyrme-type interactions are used in a 
Hartree-Fock scheme. It is found that the critical density 
for ferromagnetism decreases with temperature. On the 
other hand, in Ref. [10] the authors report no indication 
of ferromagnetic instability at any density or temperature 
based on the Brueckner-Hartree-Fock (BHF) approxima- 
tion and the Argonne V18 nucleon-nucleon (NN) interac- 
tion [ll|. The properties of spin-polarized neutron mat- 
ter at finite temperature are studied in Ref. with 
two different parameterizations of the Gogny interaction. 
The results show two qualitatively different behaviors for 
the two parameterizations. The reasons for these discrep- 
ancies must be carefully studied and their origin under- 
stood in terms of specific features of the nuclear force 
and/or the chosen many-body framework. 

Previous work on the temperature dependence of the 
EoS includes the calculations by Baldo and Ferreira who 
used the Bloch-De Doniinicis diagrammatic expansion 
[13], Brueckner-Hartree-Fock calculations with and with- 
outh three-body forces (TBF) [2|], and the predictions of 
Ref. [13] based on the Green's function method. Inves- 
tigations of the finite temperature EoS have also been 
performed within the relativistic Hartree approximation 
[l5l - l2l| . whereas the relativistic Hartree-Fock approach 
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was applied in Refs. [21j. 

The entropy per particle in symmetric nuclear matter 
has been studied in Ref. [25] within the self-consistent 
Green's function (SCGF) approach, where both particle- 
particle and hole-hole scatterings are included. The 
SCGF framework allows direct access to the single- 
particle spectral function and thus to all the one-body 
properties of the system. A most recent work by Rios 
et al. (23j addresses hot neutron matter within the same 
approch and performs a comparison with other models. 
Hot asymmetric matter and /3-stable matter have been 
studied by Moustakidis et al. [13, HI] using temperature 
and momentum dependent effective interactions. An ear- 
lier calculation with the DBHF method can be found in 
Ref. [t^. Also, Ref. [%] contains a study of isospin asym- 
metric matter based on the DBHF approach. 

The study of the many EoS-related aspects in both 
symmetric and neutron-rich matter starting from real- 
istic NN forces and within a microscopic model is still 
a considerable challenge. It is the purpose of this pa- 
per to report the first part of a comprehensive study 
of temperature dependence of nuclear and neutron-rich 
matter properties based on the Dirac-Brueckner-Hartree- 
Fock (DBHF) approach. Here, we will concentrate 
on the properties of the single-particle within the nu- 
clear medium for the following reason: although the nu- 
clear/neutron matter energy density is certainly an im- 
portant quantity, the single-particle interaction and its 
temperature dependence, which determine the one-body 
properties in the medium, are perhaps more relevant for 
non-equilibrium processes such as relativistic heavy-ion 
collisions. 

Previously, we have confronted isospin asymmetries 
[26l [27I as well as spin asymmetries effects on the 
equation of state of cold matter. In this work, we have 
extended the isospin-asymmetric matter calculation of 
Ref. [26*] to include finite temperature effects. The simul- 
taneous consideration of isospin asymmetry and temper- 
ature dependence will make our microscopic predictions 
more broadly useful and capable to reach out to the prop- 
erties of the hot environment present in the latest stage 
of a supernova collapse or in the collision of heavy nu- 
clei at intermediate energies. Starting from the present 
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baseline, we plan to address, in future work, additional 
aspects such as temperature dependence of in-medium 
effective cross sections and hot spin asymmetric matter. 

Next, after a review of the main aspects of the formal- 
ism, we will show and discuss predictions of the chemical 
potential and single-particle properties in symmetric and 
pure neutron matter. In Section IIIC we show a typical 
set of predictions for the single-neutron and single-proton 
potentials in isospin-asymmetric matter. 

Hot nuclear and neutron-rich matter are infinite 
fermionic systems where only the strong interactions 
among nucleons are taken into account. The tempera- 
tures typically considered are of a few tens of MeV, rela- 
tively small on the scale of nuclear energies. For instance, 
close to the saturation density of nuclear matter, the free 
Fermi energy, ep, is approximately 40-50 MeV, and thus 
a temperature of 20 MeV is still somewhat low, with 
T/ep ^ 0.5. 



II. FORMALISM 

Within the DBHF method, the interactions of the nu- 
cleons with the nuclear medium are expressed as self- 
energy corrections to the nucleon propagator. That is, 
the nucleons are regarded as "dressed" particles, essen- 
tially a gas of non-interacting quasi-fermions. The be- 
havior of the dressed nucleon is determined by the effec- 
tive nucleon propagator, which obeys the Dyson equa- 
tion. Relativistic effects lead to an intrinsically density- 
dependent interaction which is consistent with the contri- 
bution from TBF typically employed in non-relativistic 
approaches. The advantage of the DBHF approximation 
is the absence of phenomenological TBF to be extrap- 
olated at higher densities from their values determined 
through observables at normal density. 

In the quasi-particle approximation, the transition to 
the temperature-dependent case is introduced by replac- 
ing the zero-temperature occupation number with its 
finite-temperature counterpart, namely the Fermi-Dirac 
occupation density. In fact, it can be shown that the 
zero-temperature terms in the Bethe-Goldstone expan- 
sion where temperature is introduced in the occupation 
number only are the dominant ones fis'], thus justifying 
this simplified procedure. More specifically, one replaces 
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with the Fermi-Dirac distribution 



nFD{k,p,T) 



1 



I ^ g(e(fc,p,T)-M(p,T))/T 



(1) 



(2) 



Here T is the temperature in MeV, €{k,p,T) the single- 
particle energy, function of momentum, density, and tem- 
perature, and /i the chemical potential, to be determined. 
The angle-averaged Pauli operator is evaluated numeri- 



cally. It's given by 
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Q{q,P,P,T) = - j d{cos e){l-nFD{k,p,T))x 



x{l~npD{k\p,T)) , (3) 

for two nucleons with momenta k and fc', with relative 
and total momentum given by <f = and P — k + k', 
respectively. 

The single-particle energy, e{k,p,T), is now tempera- 
ture dependent. It can be obtained self-consistently with 
the Dirac states following the same procedure as used in 
the zero temperature case, but including Eq. (2) in the 
calculation of the single-particle potential. At each iter- 
ation of the self-consistent calculation, the normalization 
condition 
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npoik, p,T)d^k , 
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allows to extract the microscopic chemical potential, 
p{p, T). D is a degeneracy factor, equal to 4 for symmet- 
ric unpolarized nuclear matter or 2 for unpolarized neu- 
tron matter. For isospin-asymmetric matter, proton and 
neutron densities are fitted simultaneously, see Eq. (15). 

In close analogy with the T = case, the single- 
particle potential and the self-consistent nucleon-nucleon 
G-matrix are related by 

\ ^F, 1 ^ 2t+ 1 2s + 1 



I,L,S,J 



X / npD{k',p,T)G^^'^^^'^^\q{k,k'),P{k,k'))<i'k' , 
Jo 

(5) 

where /, L, S, and J are the NN system quantum num- 
bers while s and t are the single-particle spin and isospin, 
respectively. Integrations over the Fermi sea are per- 
formed using a cutoff adjusted to the density. We found 
a value of 6 fm to be sufficient for the highest densities 
and temperatures considered here. 

The single-particle energy is the sum of potential (U) 
and kinetic (K) energy contributions, 



eik,p,T) = U{k,p,T) + K{k,p,T). 



(6) 



It is parametrized using the same ansatz as in the T = 
case [29| for a spin-symmetric, rotationally invariant 
system: 



s{k,p,T) = Vk^ + {m*)^ + Uv, 



(7) 



where Uy is the time-like part of the vector potential 
and m* — m + Us, with Us the scalar potential, now 
both density and teperature dependent. This quantity 
plays an important role in temperature and momentum- 
dependent transport models of heavy-ion collisions. 
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FIG. 1: (color online) The chemical potential in symmetric 
matter as a function of the Fermi momentum at various tem- 
peratures from T = to T = 30 MeV in steps of 5 MeV. The 
chemical potential decreases with temperature. 
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FIG. 3: (color online) Effective mass in symmetric matter as 
a function of the Fermi momentum and for different temper- 
atures: T = (solid black); T = 10 (dashed blue); T = 15 
(dash-dotted green); T = 20 (dotted red). 




FIG. 2: (color online) The single-particle potential in sym- 
metric matter as a function of the momentum for three dif- 
ferent temperatures: T — (solid black); T = 10 (dashed 
green); T = 20 (dotted red). The left and right panels corre- 
spond to Fermi momenta equal to 0.9 fm~^ and 1.3 fm~^, 
respectively. 
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FIG. 4: (color online) The single-particle energy minus the 
rest mass in symmetric matter parametrized as in Eq. (7) vs. 
the momentum at different temperatures: T = (solid black) ; 
T = 10 (dashed green); T = 20 (dotted red). The left and 
right panels correspond to Fermi momenta equal to 0.9 fm~^ 
and 1.3 fm~^, respectively. 



III. RESULTS 



Once a self-consistent solution is obtained for the 
single-particle potential, and, simultaneously, for the 
chemical potential, the entropy /particle can be calcu- 
lated, along with all thermodynamic quantities. In the 
mean- field approximation, the entropy /particle is given 
by 

1 D f°° 

5 = ,o / [nFD{k,p,T) In npD{k,p,T) + 

P {^T^Y Jo 



+ (1 - npDik, p, T)) In (1 - npoik, p, T))] cPk , (8) 

that is, it has the same functional form as in a noninter- 
acting system. 



A. Symmetric Nuclear Matter 

We begin by showing the predicted chemical potential 
in symmetric nuclear matter (SNM), see Fig. 1. Each 
curve is an isotherm starting from zero temperature and 
going up to T = 30 MeV in steps of 5 MeV (/i goes down 
with increasing temperature). The effect of temperature 
is definitely larger at the lower densities and increases in 
size with increasing temperature. Our predictions are in 
fair qualitative agreement with those shown in Ref. [T^ , 
with or without the contribution of TBF, which do not 
seem to have a major effect on the chemical potential. 

Next, we show the single-particle potential in SNM, 
U{k, p, T), see Fig. 2. The effect of temperature is much 
more pronounced at low density (compare left and right 
panels) and low momenta. The general tendency is to 
turn slightly more repulsive with increasing temperature, 
although this trend becomes clear only at the higher tem- 
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FIG. 5: (color online) The entropy /particle in SNM as a func- 
tion of density and for increasing temperatures of T = 5, 10, 
15, and 20 MeV. The entropy increases with temperature. 



FIG. 6: The entropy /particle in SNM as a function of tem- 
perature at a density corresponding to a Fermi momentum of 
1.3 /m~\ 



perature. The increased repulsion is the result of a com- 
bination of effects. Temperature "smears out" the step 
function distribution, Eq. (1), so that the interaction 
probability increases (decreases) at low (high) momenta 
due to the smaller (larger) occupation probability as com- 
pared to the r = case. Although Pauli blocking is gen- 
erally reduced by temperature (which suggests increased 
attraction among the particles), the integral in Eq. (5) 
receives contributions from G-matrix elements at higher 
momenta (as compared to the zero temperature case), 
and such contributions tend to be repulsive. In the end, 
we observe a net effect that is repulsive, except at the 
lowest temperatures, and essentially negligible at high 
momenta. In Ref. J] ^^e equivalent Schroedinger opti- 
cal potentail was constructed from the Dirac self-energies 
and found to be remarkably insensitive to temperature. 

Consistent with Fig. 2, the temperature dependence of 
the effective mass is generally small, see Fig. 3, and more 
noticeable at low density. Thus, in the range of densities 
and temperatures considered here, we expect only minor 
effects on the in-medium cross sections, whose behavior 
is essentially dominated by the effective mass. A slightly 
repulsive temperature effect on the DBHF effective mass 
was found in Ref. as well. 

The single-particle energy, not including the rest mass, 
is shown in Fig. 4 at fixed density and for different tem- 
peratures. Fig. 4 confirms that the single-particle interac- 
tion is noticeably impacted only at the lowest momenta. 

We conclude this set of results by showing the en- 
tropy/particle, which is a measure of thermal disor- 
der. Entropy production in multifragmentation events in 
heavy-ion collisions is a crucial quantity in the determi- 
mation of the mass fragment distribution. The entropy 
per particle is shown in Fig. 5 as a function of density and 
for various temperatures, and in Fig. 6 it is displayed as a 
function of temperature at a fixed density (close to satu- 
ration density). The entropy increases with temperature, 
as physically reasonable, and decreases substantially with 
density. At low T, it is expected to approach a linear de- 
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FIG. 7: (color online) The chemical potential in neutron mat- 
ter as a function of the Fermi momentum at T = 0, 10, and 
20 MeV. The chemical potential decreases with temperature. 

pendence due to the fact that, for a Fermi liquid, the 
relation between S and T should be approximately 

S « ^N{T = 0)r, (9) 

in terms of the density of states at the Fermi surface. 

Although weak model dependence is not a general fea- 
ture of thermodynamic quantities, the authors of Ref. [l^ 
demonstrate that different approximations to the en- 
tropy, including the quasi-particle approximation in the 
temperature dependent BHF scheme, differ from each 
other by 10 to 20% at most. This may be due to cancela- 
tions in the difference between the single-particle energy 
and the chemical potential [2^ . 

B. Neutron Matter 

We have done a similar study of hot neutron matter 
(NM) as well. The chemical potential in NM is shown 
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namical quasiparticle" result ^2§\. 

The temperature dependence of the single-particle in- 
teraction in hot NM is comparable to, or smaller than, 
the one encountered in the SNM case. We show a repre- 
sentative case in Fig. 9. It is interesting to notice, from 
Fig. 9 and from the left panel in Fig. 2, how there is a 
tendency of the potential to become deeper at first and 
then more shallow as the temperature increases. Also, 
temperature appears to "wash out" some of the struc- 
ture in the potential. 



C. Isospin-asymmetric matter 



FIG. 8: (color online) The entropy/particle in NM as a func- 
tion of density at T = 10 (lower curve) and 20 MeV (upper 
curve) . 




FIG. 9: (color online) The single-particle potential in neutron 
matter as a function of the momentum at different temper- 
atures: T = (solid black); T = 10 (dashed blue); T = 20 
(dash-dotted green); T — 30 (dotted red). The neutron mat- 
ter Fermi momentum is equal to 0.9 fm~^. 



in Fig. 7 as a function of the Fermi momentum. (Of 
course, the same Fermi momentum corresponds to half 
the density as compared to SNM.) The microscopic chem- 
ical potential, as obtained from the normalization condi- 
tion Eq. (4), was shown to be in good agreement with 
the macroscopic one, obtained from the bulk properties 
through the derivative of the free energy density ^2^] . 

The entropy in NM is displayed in Fig. 8. As for the 
case of SNM, discrepancies between predictions from dif- 
ferent models have been found to be small, especially 
those arising from the use of different NN potential mod- 
els ^2§\. In fact, our predictions (based on the Bonn 
B potential) are close to those shown in Ref. [23] either 
with CDBONN or Argonne V18 [HI. Concerning dif- 
ferent many-body approaches, there are indications that 
contributions to the entropy from dynamical correlations 
(which would fragment the quasiparticle peak) are small. 
Hence, different predictions tend to approach the "dy- 



The isospin splitting of the nucleon mean field in 
isospin-asymmetric matter is a topic of great contem- 
porary interest. It gives rise to the symmetry potential, 
a quantity which plays a crricial role in simulations of 
heavy-ion collisions (see Ref. (2^ and references therein) . 

We have incorporate temperature dependence in 
our calculation of isospin-asymmetric nuclear matter 
(lANM) [2^, [23] ■ Asymmetric nuclear matter can be 
characterized by the neutron density, pn, and the pro- 
ton density, pp. In infinite matter, they are obtained by 
summing the neutron or proton states per volume (up to 
their respective Fermi momenta, kp or kp) and applying 
the appropriate degeneracy factor. The result is 



Pi 



37r2 



(10) 



with i = n or p. It is more convenient to refer to the 
total density p — Pn+ Pp and the asymmetry (or neutron 
excess) parameter a = BjlzEi.^ Clearly, a=0 corresponds 
to symmetric matter and a=l to neutron matter. In 
terms of a and the average Fermi momentum, kp^ related 
to the total density in the usual way. 



P - 



3^' 



(11) 



the neutron and proton Fermi momenta can be expressed 
as 



k-^^kpil + a) 



1/3 



and 



kp = kpi^ — a) 



1/3 



(12) 



(13) 



respectively. The formulas given in Section II are ex- 
tended to include two types of fermions. For instance. 
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)r,r'{q,P,P,a,T) ^ 2 J 



d{cos 6){l—npp,^T{ki,p, a,T))x 



X (1 - nFD,T'{k2,p,ai,T)) , 



(14) 
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FIG. 10: (color online) Single-neutron (upper panel) and 
single-proton (lower panel) potentials vs. momentum (in units 
of the Fermi momentum, which is equal to 1.4/m~^). In 
the upper panel, the three higher curves are the predictions 
for at T=0 (short dash), T=10 MeV (long dash), and 

T=20 MeV (solid). The three lower curves are the predic- 
tions in isospin symmetric matter at the same temperatures. 
The predictions become more repulsive with increasing tem- 
perature, as it can be best seen at the lowest momenta. In 
the lower panel, the predictions for the a— 0.3 case are below 
those corresponding to a—0. 



where r, t' =n or p. The normaUzation condition, for each 
type of nucleon, becomes: 

1 /"^ 

Pt = 2 / npD,T{k,p,a,T)d^k . (15) 

(.^""J Jo 




FIG. 11: (color online) Similar to Fig. 10, but for a=0.5. 
Similarly, the single-proton or neutron potentials become 

Urik,p,a,T)=Y,Urr' , (16) 
t' 

where each Utt' piece refers to a particular type of nu- 
cleon, and thus looks like Eq. (5) without the summation 
over isospin states. The calculation then proceeds to de- 
termine, self-consistently, the parameters of both protons 
and neutrons. 

In Figs. 10-11, we show a representative set of re- 
sults for the single-neutron and proton potentials. In 
both figures, the total density corresponds to an aver- 
age Fermi momentum equal to 1.4/m~^. The asym- 
metry parameter is equal to 0.3 and 0.5 in Fig. 10 and 
Fig. 11, respectively. The upper (lower) panel displays 
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the neutron (proton) potential. The potentials are dis- 
played as a function of the momentum (in units of the 
Fermi momentum) for three different temperatures. For 
the neutron case, the three lower curves correspond to 
a—0, whereas the three upper ones are the predictions 
for a~0.3. In both sets, the potential becomes slightly 
more repulsive with increasing temperature, something 
which is noticeable only at the lowest momenta, as ob- 
served previously. For the proton case, similar consid- 
erations apply concerning temperature effects, but here 
the a=0 potentials are above the a— 0.3 predictions. Fig- 
ure 11 displays a similar scenario, but for a=0.5. Thus, 
Figs. 10-11 show the effect of isospin asymmetry at each 
temperature. As expected [13], the neutron and the pro- 
ton potentials "move away" from the symmetric matter 
predictions in opposite directions. The effect of isospin 
asymmetry is larger than the one of temperature. 

IV. SUMMARY AND CONCLUSIONS 

We used the Dirac-Brueckner-Hartree-Fock method 
extended to finite temperatures to predict single-particle 
properties in hot SNM and lANM. For temperatures up 
to a few tens of MeV the effect of temperature is small 
except at low densities and momenta. Due to suppres- 
sion of Pauli blocking, which is most important at low 
momentum, some temperature-induced modification of 
(low-energy) in-medium cross sections could be expected. 



This is an interesting point we will explore in future work. 
We also discussed the entropy/particle. When compared 
with other predictions in the literature, our results con- 
firm the weak model dependence of this quantity. 

As usual, we adopt the microscopic approach for our 
nuclear matter calculations. Concerning our many-body 
method, we find DBHF to be a good starting point to 
look beyond the ground state of nuclear matter, which it 
describes successfully. The main strength of this method 
is its inherent ability to effectively incorporate crucial 
TBF contributions [23| yet avoiding the possibility of 
inconsistency between the parameters of the two- and 
three-body systems. 

We have focussed on the one-body properties of the 
system, which are of great relevance for the study of en- 
ergetic heavy-ion collision dynamics. The extension of 
our isospin-asymmmetric nuclear matter calculation, ini- 
tiated in Ref. ^2d\ and finalized in Ref. jl^], to include 
the effect of finite temperature will enable us to consider 
additional aspects such as, for instance, proton fraction 
in hot beta-stable matter. 
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